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ABSTRACT 

Medium effects which include color screening and Landau-Pomeranchuk-Migdal 
suppression of induced radiation are discussed in connection with the equilibration of 
dense partonic system in ultrarelativistic heavy ion collisions. Taking into account 
of these medium effects, the equilibration rate for a gluonic gas are derived and the 
consequences are discussed. 

1. Introduction 

Strong interactions involved in hadronic collisions can be generally divided into 
two categories depending on the scale of momentum transfer Q 2 of the processes. 
When Q 2 ~ Aq CD , the collisions are nonperturbative in QCD and are considered soft. 
These interactions can be approximated by some effective theories in which partons 
inside a nucleon cannot be resolved and nucleons interact coherently by the exchange 
of mesons or soft Pomerons. These kind of coherent interactions will result in the 
collective excitations as observed in experiments at low and intermediate energies, 
y/s < a few GeV. On the other hand, if Q 2 is much larger than Aq CD , parton model 
becomes relevant and they interact approximately incoherently. Due to the high 
Q 2 , parton interactions can be calculated via perturbative QCD (pQCD), given the 
initial parton distribution functions inside a nucleus. At high collider energies as 
RHIC and LHC, the hard processes of parton interactions become dominant and 
heavy ion collisions are generally very violent!. In these violent collisions, it is quite 
likely that colors are liberated during the collisions and we can start with a color 
deconfmed system. 

In the last few years, many efforts have been made to estimate the initial parton 
density by perturbative calculations!®!. Though the results depend on parton struc- 
ture functions inside a nucleus and the modeling of the underlying soft interactions, 
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a general consensus is that a dense partonic system will be produced during the early 
stage of ultrarelativistic heavy ion collisions. However, the system is not necessarily 
in both thermal and chemical equilibrium. The initial system is totally dominated by 
gluons and the quark and anti-quark densities are far below their chemical equilibrium 
values. Numerical simulations of parton cascade can be made with certain simpli- 
fications. However, the complexity of these calculations makes it difficult to obtain 
a clear understanding of the different time scales on various parameters and model 
assumptions. Most importantly, medium effects and interferences are difficult to in- 
corporate. In this talk, I would like to describe two medium effects. Color screening 
is caused by the interaction of a color charge with the medium which gives a natu- 
ral infrared cut-off to regularize the parton cross sections inside the dense system!. 
Another medium effect is the so-called Landau-Pomeranchuk-Migdal effect of the re- 
duced radiation by multiple scattering through the medium!. These medium effects 
can be utilized to obtain a parameter-free set of equations, based on perturbative 
QCD in a dense partonic medium, that describes the evolution of quark and gluon 
distributions towards equilibrium!. 



2. Color Screening 

Let us first look at the color screening in the initially produced parton system. 
Following the standard calculation(in Coulomb gauge) of screening in the time-like 
gluon propagator in a medium of gluonic excitation, we have the screening massi, 

& = lirn f d 3 k^- q ■ V k /(k). (1) 

7T |q|->0 J q • k 

Now, instead of using the Bose-Einstein distribution for the thermalized case, we 
simply relate the phase space density /(k) to the initial gluon distribution calculated 
from pQCD. At high energies we find that the transverse and longitudinal screening 
lengths are very close. We assume that the gluon momentum distribution is similar 
to a thermal distribution 

/„(*) = \e~ u - kt \ (2) 

with X g and T characterizing the thermalization of the gluon gas. With this distri- 
bution, one obtain the effective color screening mass 

A = X 9 g 2 T 2 . (3) 

We will use this effective color screening mass as an infrared cut-off for parton inter- 
action cross sections inside a non-equilibrium gluonic gas. 



3. Induced Radiation 

The leading processes which contribute to the gluon equilibration are those with 
induced gluon radiation. Let us consider the simplest case of induced radiation from a 
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two-quark scattering. The Born amplitude for two-quark scattering (pi, kj) — > (jpf, kf) 
through one gluon exchange is, 



u(pf)ipu(pi)u{kf)Yu{ki) 



M el = ig'T^T^, fe )2 . ( 4 ) 

where A, A', B, and B' are the initial and final color indices of the beam and target 
partons. The corresponding elastic cross section is, 
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where C^' = Cp/2N = 2/9 is the color factor for a single elastic quark-quark scat- 
tering and s,u, and t are the Mandelstam variables. 

Taking into account the dominant contribution, the radiation amplitude from a 
single scattering is, 

a a _ ■ Mel „ 

M rad = * a T a U ^ 
I AA' 1 BB' 



Tlx ~ 2ige L ■ 



k\ + (cu " 



n A ,[T a ,T b } BB ,, (6) 



where M. e \ is the elastic amplitude as given in Eq. ([|), and 72. i is defined as the 
radiation amplitude induced by a single scattering. For later convenience, all the 
color matrices are included in the definition of the radiation amplitude 72i. With the 
above approximations, we then recover the differential cross section for induced gluon 
bremsstrahlung by a single collision as originally derived by Gunion and Bertsch i, 
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dtdyd 2 k± dt dyd 2 k± 
where the spectrum for the radiated gluon is, 

dn^ _ 1 — -2 _ C A a s q\ 
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(7) 



dyd 2 k ± 2(27r) 3 C'2 ) tt 2 ki(q ± -k ± ) 2 ' 

In the square modulus of the radiation amplitude, an average and a sum over initial 
and final color indices and polarization are understood. The above formula is also 
approximately valid for induced radiation off a gluon line, except that the color factor 
Cj^ in the elastic cross section has to be replaced by 1/2 for gq and 9/8 for gg 
scatterings. 

One nonabelian feature in the induced gluon radiation amplitude, Eq. (|6|), is the 
singularity at = q_i_ due to induced radiation along the direction of the exchanged 
gluon. For kj_ we note that the induced radiation from a three gluon vertex 

can be neglected as compared to the leading contribution l/k\. However, at large 
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k± ^> q±, this three gluon amplitude is important to change the gluon spectrum to a 
l/k\_ behavior, leading to a finite average transverse momentum. Therefore, q± may 
serve as a cut-off for kj_ when one neglects the amplitude with the three gluon vertices 
as we will do when we consider induced radiation by multiple scatterings in the next 
section. 



4. EFFECTIVE FORMATION TIME 

The radiation amplitude induced by multiple scatterings has been discussed in 
Ref . |7j . We here only briefly discuss the case associated with double scatterings. We 
consider two static potentials separated by a distance L which is much larger than 
the interaction length, For convenience of discussion we neglect the color indices 
in the case of an abelian interaction first. The radiation amplitude associated with 
double scatterings is, 



7^ ED = ze 
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where p = {p°f,p z ,P±) is the four-momentum of the intermediate parton line which 
is put on mass shell by the pole in one of the parton propagators, X\ = (0,xi), and 
X2 = (^2,^2) are the four-coordinates of the two potentials with t 2 = {z-2 — z-\)/v = 
Lp°/p z . We notice that the amplitude has two distinguished contributions from each 
scattering. Especially, the diagram with a gluon radiated from the intermediate line 
between the two scatterings contributes both as the final state radiation for the first 
scattering and the initial state radiation for the second scattering. The relative phase 
factor k ■ (x 2 — x\) = uo{l/v — cos 6)L then will determine the interference between 
radiations from the two scatterings. If we define the formation time as 



2u 

uj(1/v — cos9) kj 
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then Bethe-Heitler limit is reached when L ^> r(k). In this limit, the intensity of 
induced radiation is simply additive in the number of scatterings. However, when 
L ^ r(k), the final state radiation from the first scattering completely cancels the 
initial state radiation from the second scattering. The radiation pattern looks as if the 
parton has only suffered a single scattering. This is often referred to as the Landau- 
Pomeranchuk-Migdal (LPM) effectS The corresponding limit is usually called fac- 
torization limit. 

The extrapolation to the general case of m number of scatterings gives us the 
radiation amplitude in QCD, 

? k m 
K m = i2g—±T£ A , ■ ■ ■ TT mA , m E {T am • • • [T a \T b ] • • • 7*)^, e** (11) 
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which contains m terms each having a common momentum dependence in the high 
energy limit, but with different color and phase factors. The above expression should 
also be valid for a gluon beam jet, with the corresponding color matrices replaced by 
those of an adjoint representation. The spectrum of soft bremsstrahlung associated 
with multiple scatterings in a color neutral ensemble is, similar to Eq. (H), 



dn^ 
dyd 2 k± 
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dyd 2 kj 



(12) 



where = (CV/2iV) m is the color factor for the elastic scattering cross section 

without radiation. C m (k), defined as the "radiation formation factor" to characterize 
the interference pattern due to multiple scatterings, can be expressed as, 
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(13) 
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where the color coefficients are defined as 

(J. . — ^pa m _ _ _ ]rpb rpcii^ _ _ _ rpairpai _ _ _ ^paj J^^j . . . l J' a m.j (14) 

If we average over the interaction points x« according to a linear kinetic theory, 
we find that an effective formation time in QCD can be defined as, 

Ca 2 cosh y 



tqcd(&) = r 2 r(k) 



2C 2 k 



(15) 



which depends on the color representation of the jet parton. The induced gluon 
radiation due to multiple scattering will be suppressed when the mean free path A is 
much smaller than the effective formation time. 



5. Equilibration Rates 

The most important reactions for establishing gluon equilibrium are gg <-> ggg. 
Elastic scattering processes, on the other hand, are crucial for maintaining local ther- 
mal equilibrium. Multi-gluon radiation is presumably suppressed by color screening, 
while radiative processes involving quarks have smaller cross sections due to QCD 
color factors. 

The evolution of the gluon density n according to reactions mentioned above can 
be described by a rate equation. Adding the equation for energy conservation assum- 
ing only longitudinal expansion we end up with a closed set of equations determining 
the temperature T(r) and the gluon "fugacity" \ g {j) = n/n eq (T) as a function of the 
proper time r, 

R 3 (l-X g ), (16) 
const. (17) 
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This set of evolution equations is completely controlled by the gluon production rate, 
R3 = ^ 3 n. 

Taking into account of LPM effect, the cross section for gg — > ggg processes can 
be written as, 

<PqX 3 d?k ± = ^^^ k±C ° Shy9{X ~~ T ^ {k))d{E ~ k ^° Shy) > (18) 

where tqcd(^) is given by Eq. (fT5P, the second ^-function is for energy conservation. 
The gluon density distribution induced by a single scattering is given by Eq. (H) which 
must be regularized by the color screening mass //p in Eq. (|3|). Using the elastic cross 
section of gluon scattering regularized also by //£>, we obtain a fugacity independent 
mean free path 

A" 1 = a 2 ^ 2 n = §aia,T, (19) 

Using these values we evaluate the chemical gluon equilibration rate R3 = \na 3) 
as defined in Eq. (0), numerically. This rate scales with the temperature linearly 
but is a complicated function of the gluon fugacity. The result can approximated by 
an analytical fit, 

J R 3 = 2.1a s 2 T(2A 9 -A2) 1/2 , (20) 

which will be used in solving the time dependent rate equations. The rate we thus 
obtained has a nonlinear dependence on the fugacity due to the inclusion of the LPM 
effect and the effective color screening mass. If LMP effect were not included, the 
rate whould be much larger, thus would lead to a fast equilibration. 
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